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I. INTRODUCTION 



A number of recent cosmological observations, e.g., Type la Supernovae baryon acous 



tic oscil 
diation 



ations [2J, large scale structure (LSS) |3|, cosmic microwave background (CMB) ra- 

4- 



6], and weak lensing suggest the accelerated expansion of the current universe. 
To explain the late time cosmic acceleration, there exist two main approaches: One is the 
introduction of the so-called "dark energy" (for reviews on dark energy, see, e.g., fl-Q) and 



the other is the modification o 
Recently, "teleparallelism" 



gravity such as f(R) gravity [15|-|23|. 



24J has attracted much attention as it can be considered as 



an alternative gravity theory to general relativity. Teleparallelism is formulated with the 
Weitzenbock connection, so that its action consists of the torsion scalar T instead of the 
scalar curvature R in general relativity with the Levi-Civita connection. It has been shown 
that by introducing a scalar with the non- minimal coupling to gravity in teleparallelism 25| . 
the late time cosmic acceleration can be achieved. Moreover, similar to f(R) gravity, the non- 
linear generalization of the torsion scalar T, i.e., f(T) gravity, can account for inflation 26] 



in the early universe as well as the cosmic acceleration in the late time 



27, 



281. 



Various 

aspects on f(T) gravity have been widely investigated in the literature (see, e.g., [ijj] and 
the references therein). 

On the other hand, according to astrophysical observations, it is well known that there 
exist magnetic fields with the strength ~ G and the coherence scale 1-10 kpc. Also 
in clusters of galaxies, large-scale magnetic fields are observed, whose strengths are 1CT 7 - 
10 -6 G and the coherence scales are estimated as 10 kpc-1 Mpc. However, the origins of these 
cosmic magnetic fields, in particular the large-scale magnetic fields in clusters of g alaxies have 
not been well understood yet (for reviews on cosmic magnetic fields, see, e.g., 29]). There are 
several generation mechanisms of the cosmic magnetic fields, such as those from astrophysical 



processes based on the plasma instability 



30 



31], cosmological phase transitions 32], and 



matter density perturbations before or at the recombination epoch 33]. Indeed, it is not so 
easy for these mechanisms to generate the large-scale magnetic fields observed in clusters 
of galaxies only with the adiabatic compression and without any secondary amplification 
mechanism as the galactic dynamo 34|. Thus, the most natural mechanism to produce the 
large-scale magnetic fields is considered to be electromagnetic quantum fluctuations during 
inflation 35] , because the scale of the electromagnetic quantum fluctuations can be extended 
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to that larger than the Hubble horizon by inflation. 

In Quantum Electrodynamics (QED) in the curved space-time, there can appear a non- 
minimal coupling of the scalar curvature to the electromagnetic field owing to one-loop 
vacuum-polarization effects |J 



so that the conformal invariance of the electromagnetic 
field can be broken by this coupling. This can yield the quantum fluctuation of the electro- 
magnetic field during inflation, resulting in the large-scale magnetic field at the present uni- 
verse HQ. Such a breaking mechanism of the conformal invariance of the electromagnetic 
field is necessary to generate the quantum fluctuation of the electromagnetic field, because 
the Maxwell theory is conformally invariant and the Friedmann-Lemaitre- Robertson- Walker 
(FLRW) space-time is conformally flat 38I] 1 . Consequently, a lot of breaking mechanisms of 
the conformal invariance of the electromagnetic field have been explored (for a list of these 



4J-|52|). 



breaking mechanisms, see, e.g., reviews in [29| and references in 

In this paper, motivated by both astrophysical and cosmological observations, we study 
the generation of large-scale magnetic fields from inflation in teleparallelism. In particular, 
we introduce a non-minimal gravitational coupling of the torsion scalar T to the electro- 
magnetic field by analogy with such an interaction between gravity and electromagnetism 
in general relativity. As an illustration, we demonstrate that for the form of the coupling to 
be a power-law type, the magnetic field with its current strength of ~ 10 _9 G on IMpc scale 
can be generated. 



It should be remarked that for example, in Ref. 44] Ratra has investigated the case 
that the gauge kinetic term is coupled to the inflaton field. In this work, however, the 
observation that there can be a spectator field evolving during inflation is used essentially. 
This possibility has been scrutinized in a number of different works, e.g., Refs. 53M55] by 



Giovannini. In particular, in Ref. [53( a scale-invariant spectrum during the conventional 
inflation has been demonstrated in a specific model where the gauge coupling is not a function 
of the inflaton (in the latter case the flatness of the potential might be spoiled). We use 



units of &b = c = H — 1 and denote the gravitational constant 8irG by n 2 
the Planck mass of M Pl = G~ 1/2 = 1.2 x 10 19 GeV. 



/Mpi 2 with 



1 It should be noted that this is true only for the flat FLRW space-time, but not for the FLRW background 
with spatial curvature, e.g., an open FLRW universe 1391 . In addition, the breaking of the conformal 
flatness during inflation has also been studied in Ref. [40(. Furthermore, there exist arguments in terms 
of the back reaction effect of the magnetic field generated during inflation 4ll443[. 
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The paper is organized as follows. In Sec. II, we explain the fundamental formulations in 
teleparallelism. In Sec. Ill, in a non-minimal /(T)-Maxwell theory, where I{T) is an arbitrary 
function of the torsion scalar T, we investigate the generation of large-scale magnetic fields 
in inflationary cosmology. In Sec. IV, for a concrete model of a power-law type coupling 
between the torsion scalar and the Maxwell field, we estimate the current strength of the 
large-scale magnetic field. Finally, conclusions are presented in Sec. V. 



II. TELEPARALLELISM 



We adopt orthonormal tetrad components eA{x^) in teleparallelism, where an index 
A runs over 0,1,2,3 for the tangent space at each point x M of the manifold. The re- 
lations between the metric and orthonormal tetrad components are given by g^ v = 
VABefev, where /i and v are coordinate indices on the manifold and run over 0,1,2,3. 
Hence, e A form the tangent vector of the manifold. We define the torsion and con- 
torsion tensors as = e p A (d^ - d v ef) and p = - (1/2) (T» u p - p - T p , 
respectively. Using these tensors, we construct the torsion scalar T = S p ^ u T p pu with 
S f / U = (1/2) (K^ p + 5% T a \ - 8 v p T Q ^ Q ). In general relativity, the Einstein- Hilbert action 
consists of the scalar curvature R. However, in teleparallelism the torsion scalar T is used 
to represent the teleparallel Lagrangian density. As a result, the action in teleparallelism is 
described by 

S Tel = J d 4 x\e\ [T/ (2k 2 ) + £ M ] , (2.1) 

where |e| = det (e^) = y/—g and Cm is the Lagrangian of matter. The variation of the action 
SVci with respect to the vierbein vector fields e\ leads to the gravitational field equation 2^], 
given by (1/e) 0„ (eS/ v ) - e x A T p pX S p ^ + (1/4) e A T = (k 2 /2) e p A T^ p \ where T^ M \ V is the 
energy-momentum tensor of matter. 

We take the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, whose metric 
is given by ds 2 = dt 2 — a 2 (t)dx 2 = a 2 (r]) (—drj 2 + dx 2 ) with a the scale factor and rj the 
conformal time. In this space-time, g^ u = diag(l, —a 2 , —a 2 , —a 2 ) and the tetrad components 
become e p = (l,a, a, a). With these relations, we find that the exact value of the torsion 
scalar is described by T = —QH 2 , where H = a/ a is the Hubble parameter with the dot 
being the time derivative of d/dt. 
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III. NON-MINIMAL /(T)-MAXWELL THEORY 



In this section, we consider a non-minimal /(T)-Maxwell theory and examine the gener- 
ation of large-scale magnetic fields in inflationary cosmology. 



A. Model of the electromagnetic sector 



The action describing a non-minimal /(T)-Maxwell theory is given by 



S= I d*x\e\{-h(T)F^) , (3.1) 



where I{T) is an arbitrary function of the torsion scalar T and F^ u = d^A u — d u A^ with 
the U(l) gauge field is the electromagnetic field-strength tensor. It follows from the action 
in Eq. (13. ip that the electromagnetic field equation is derived as 

■ -.d,[^I(T)F^] =0. (3.2) 



In the FLRW background with the Coulomb gauge A (t, x) = and &*Aj(t, x) = 0, the 
equation of motion for the U(l) gauge field is written as 

A\(t, x)+(h+ j Mt, x) - 1 ( A A { (t, x) = 0, (3.3) 

(3) 

where A is the Laplacian in three-dimensional space. 

We remark that although the electromagnetic field coupled (preferentially) to the axial 
vector part of the torsion tensor is more natural, we do not consider it in this study since the 
resulting large-scale magnetic field is too small. In addition, we emphasize that the strength 
of the generated magnetic field is more related to the nature of the coupling / to the gauge 
kinetic term than to the torsion scalar itself. 

It is important to explicitly state that the key element to generate the large-scale mag- 
netic fields in the quasi-de Sitter phase of expansion, i.e., inflation in the early universe, is 
the scalar coupling between / and the kinetic term of the electromagnetic field, as seen in 
Eq. (13. ID . In general, the coupling term I(r]) may be a function of various scalar degrees of 
freedom existing in a model, e.g., the inflaton or the dilaton field or a dynamic gauge cou- 
pling. As a result, I(r)) can be a function of a spectator field evolving during the inflationary 
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epoch. In this case, there is no connection between the evolution of / and the gauge cou- 
pling. Therefore, the physical features of the various models are different. In other words, in 
principle / can be an arbitrary function of some non-trivial background fields. On the other 
hand, in the case of bouncing models some of these ideas are preferentially realized, whereas 
other models are consistent with the standard inflationary paradigm. In this work, instead 
of concentrating on a specific mechanism for inflation, we execute a model-independent anal- 
ysis on the generation of large-scale magnetic fields through the breaking of the conformal 
invariance of the electromagnetic field due to the coupling of / in in — (1/4) IF 11V F^ U in 
Eq. flHH). 

We also describe the realization of inflation. In this work, we suppose that the generic 
slow-roll inflation is realized without identifying the specific mechanism to lead to inflation. 
There are various possibilities to realize inflation. For example, one can introduce the inflaton 



coupling to the electromagnetic field |44l. |48| . which can be considered to be the dilaton 46], 



or both inflaton and dilaton fields |47j. Namely, if the dilaton field is not responsible for 
inflation, the coupling / can be a function of the inflaton field. Furthermore, in these cases 
forms of the inflaton potential are assumed to be flat enough to realize the slow-roll inflation, 
namely, the quasi exponential expansion of the universe. Concrete demonstrations have been 



investigated in Refs. [53l-l55|. 



B. Quantization 



We now execute the quantization of A^(t,x). From the action of the electromagnetic 
fields in Eq. (13. ip . we find that the canonical momenta conjugate to A^it, x) become tiq = 0, 
7Tj = Ia(t)Ai(t, x). The canonical commutation relation between Ai(t, x) and 7ij(t,x) is 
given by [ Ai(t, x), TTj(t, y) ) = if d 3 k(2Tr)~ 3 e lk ^~ v ^ (5ij — kikj/k 2 ) with k being the co- 
moving wave number and k = \k\. By imposing this relation, Ai(t,x) is described as 
Ai{t,x) = Jrf 3 A:(27r)- 3 / 2 ^ (7=12 \b{k,a)e l {k,a)A{t,k)e ik -- + b^k,a)e*(k,a)A*{t,k)e- ih ^], 
where €i(k,a) (cr = 1,2) stand for the two orthonormal transverse polarization vec- 
tors, and U^(k,a) is the annihilation (creation) operator, satisfying the relations 



b{k,a),ti{k',a') 



0. 



= 5 a ^5 3 (k - k') and b(k, a),b(k' , a' 
It follows from Eq. fl3T3|) that the Fourier mode A(k, t) obeys A(k, t) + [H + I/I j A(k, t) + 
(fc 2 / a 2 ) A(k , t) = together with the normalization condition, A(k,t)A*(k,t) — 
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A(k, t)A*(k, t) = ij (la). If we use the conformal time rj, this equation is rewritten as 
A"(k,7]) + (r/I)A'(k,rj) + k 2 A(k,7]) = 0, where the prime denotes the derivative in terms 
of 77 as d/drj. 



C. Procedure to obtain analytic solutions 



With the WKB approximation on subhorizon scales and the long wavelength approxi- 
mation on superhorizon scales and matching these solutions at the horizon crossing, it is 
possible to acquire an analytic solution for this equation 49] approximately. In this case 
of the exact de Sitter background, we find a = 1/ (—Hrj) with H being the Hubble param- 
eter during the de Sitter expansion, and —kr] = 1 at the horizon-crossing when H = k/a. 
For subhorizon (superhorizon) scales, we have k\r]\ ^> 1 (k\i]\ <C 1). This is considered to 
be sufficiently well defined also for the general slow- roll inflation, i.e., nearly exponential 
inflation. 

Provided that in the short-wavelength limit of kj (aH) ^> 1 the vacuum asymptoti- 
cally approaches the Minkowski vacuum, the WKB subhorizon solution reads Ai n (k,rj) = 
(l/V2k^) J-VZe-^. 

While, with the long- wavelength expansion in terms of k 2 , we can have the solution on 
superhorizon scales A out (k,r)). By matching this solution with the above WKB subhorizon 
solution at the horizon crossing time rj = 77^ ~ 1/k, we obtain the lowest order approximate 
solution of A OVLt (k,rj) [49| as 



A out {k,r]) = A 1 {k) + A 2 {k) 
1 



IB 



I(fj) 



A x {k) = 

A 2 (k) EE 



j-1/2 



2k 
1 



I 2 J + lkI 





Q -ikri 











2k \2 



1 



ikl 



-ikrj 



(3.4) 
(3.5) 
(3.6) 



V=Vk 



We neglect the decaying mode solution, which is the second term of the right-hand side of 
Eq. (13. 4p . Equations (13. 5 j) and (13. 6p lead to \A(k,r))\ 2 at the late times, given by 

2 



\A(k, V )\ 2 = \A 1 (k)f 



2kI( Vk 



1 I'(Vk 



'Ik 



(3.7) 



2kI{ Vk ) 

Here, we have supposed the instantaneous reheating after inflation and therefore, ?7r is 
considered to be the conformal time at the reheating stage. By using the comoving mag- 
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netic field Bi(t,x), the proper magnetic field is expressed as /?j proper (t, x) = a~ l Bi(t,x) = 
a~ 2 £ijkdjAk(t,x), where e^* is the totally antisymmetric tensor with €123 = 1. Accord- 
ingly, we find that the spectrum of the magnetic field is described as \B proper (k, 77) | 2 = 
2 (fc 2 /a 4 ) \A(k,r])\ 2 = 2 (/c 2 /a 4 ) |t4i(/c)| 2 , where we have taken into account the factor 2 origi- 
nating from the two degrees of freedom for the polarization. In the Fourier space, the energy 
density of the magnetic field becomes pB(k,rj) = (l/2)|/? proper (/c, r])\ 2 I(r]). With multiplying 
this by the phase-space density 47r/c 3 /(27r) 3 , we derive the energy density of the generated 
magnetic field per unit logarithmic interval of k as 

PbM S l^\ B ^ k ^)Mv) = (3 ' 8) 
Consequently, the density parameter of the magnetic field per unit logarithmic interval of k 



and its spectral index are given by 



42| 



fijj(M) = — — — — - = -j-j- — - — - — I (77 , 3.9 

^ _ d\nQ B (k) _ 1 | dhxk^jk)] 2 
d In k dink 

respectively, where, p 7 (?7r) = N c s (7r 2 /30) Tr js6| is the energy density of radiation at the 
reheating stage with the reheating temperature Tr, <2r is the scale factor at 77 = t/r and N e g 
is the effective massless degrees of freedom (e.g., for photons, 2) thermalized at the reheating 
stage. 

IV. LARGE-SCALE MAGNETIC FIELD GENERATED IN TELEPARALLELISM 
A. Current strength of the magnetic field 

For the purpose of analyzing the strength of the magnetic field quantitatively, we examine 



the case of the specific form of 1(f)), given by 49] 



where 77* is some fiducial time at the inflationary stage, 7* is the value of I at 77 = 77*, and j3(> 
0) is a positive constant, whose positivity makes I increase monotonically during inflation. 
For this form of J, l^l 2 in Eq. (E3D reads A;|^4.! | 2 = (1/2/(77*)) \1 - (P + 2i) / [2(0 + 1)]| 2 = 
Aj (2/(77*)), where A{= 0(1)) is a constant of the order of unity. By plugging this relation 
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into Eqs. ( 13. 9 j) and ( I3.10p . the density parameter of the magnetic field at the present time rj 
is expressed as ^b(A;, r/o) = [A; 4 / (T R a R )] [15.4./ (2A^" e ff 7T 4 /*)] (?7fc/?7*)^ with = 4 — /?, where 
we have used oc k 13 and 1(770) = 1- If ^* is very small and the spectrum is nearly 

scale-invariant, i.e., ~ 4, the resultant amplitude of the large-scale magnetic field becomes 
large. From Eq. ( I3.9p . we find that the current density parameter of the magnetic field is 
described by 49] 



N, 



cff 



R 



-^r) 



4-/3 



(4.2) 



1080 \M Pl , 

Here, we have used the relation a R ?7 R ~ i? R 2 with a R and if R being the scale factor and 
the Hubble parameter at the reheating stage, respectively, and the Friedmann equation 
3if R / = p 7 (^r)/M| 1 at the reheating stage, where Mp\ = M P i/\/8tt = 1/k. We can further 
rewrite the term (— krjn) as (so | 



k 



R 




5.1 x 10 



(4.3) 
(4.4) 



In deriving Eq. ( 14. 3p . we have used H 1 = 3.0 x 10 3 h 1 Mpc and T oc a 1 , which leads 
to (oq/or) — (Tr/Tq). Moreover, in analyzing Eq. ( 14. 4p . we have adopted the Friedmann 



equation 3H^/ = p 7 (r7 R )/M 



1(T 29 K [56J with h 
is given by |£>(?7o)| 2 



wit 



h Av(//r) = iVeff (vr 2 /30) T 4 , T = 2.73 K and H = 2.47/t x 



0.7 



57 



Since the current amplitude of the magnetic field 
2fi s (?7o, fc) p 7 (r/ ), with p 7 (r/ ) — 2 x 10 _51 GeV 4 and 



1 G = 1.95 x 1(T 20 GeV 2 we find |50j 



\B( Vo ,L)\ = 2.7 



7.2 



1/8/8 



(5.1)% 



x ^-56+51/3/4^-4)78 J^jM 



Mi 



pi 



/3/4 



[Mpc] 



0/2-2 



G. 

(4.5) 



We note that the reheating temperature T R is described by using the Hubble pa- 
rameter at the end of inflation, namely, instantaneous reheating stage, iJ R as T R = 
[90/ (87r 3 A?efr)] 1//4 ^/M-piH-R. Furthermore, there exists the upper limit of iJ R from tensor 
perturbations. With the Wilkinson Microwave Anisotropy Probe (WMAP) five year data 
in terms of the anisotropy of the CMB radiation j^], we have i?R < 6.0 x 10 14 GeV 59]. 
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B. Estimation of the current strength of the large-scale magnetic field 



We suppose that power-law inflation occurs, in which the scale factor is given by 

a = a (^-} , (4.6) 



to, 

with p ^> 1, where ao and to are constants. The larger the value of p is, the closer power-law 
inflation goes to exponential inflation. In this case, with the relation rj — f (1/a) dt, we get 

f = [a t (p-l)(- V )r 1 ^ . (4.7) 
to 

We examine the case of a power-law type coupling as 

KT) =(!)", (4.8) 

where T is a current value of T and 0) is a non-zero constant. In this case, by using 
T = —QH 2 , H = p/t and Eq. gZD, we obtain 

7(77) = (-6/T0)" (p/t ) 2n Mo (p - l)f n/ ^ l) (- V ) 2n ^ . 

By comparing this equation with Eq. (14. li . we acquire 

In 

= -^-. (4.9) 
p — 1 

For iV cff = 100, H K = 1.0 x 10 14 GeV (T R = 8.6 x 10 15 GeV), L = IMpc, A=l, I(vr) = Hvo), 
and /3 = 4.2, which can be realized for p = 10 and n = —18.9, we have 

\B(rj ,L = lMpc)| = 2.5 x 10^ 9 G . (4.10) 

Similarly, for the above values except H R = 1.0 x 10 10 GeV (T R = 8.6 x 10 13 GeV) and 
/3 = 4.6, met for p = 10 and n = —19.7, we obtain 

\B(r] ,L = lMpc)| = 2.3 x 10~ 9 G . (4.11) 

Here, it should be mentioned that in order to demonstrate the estimation of the generated 
magnetic field strength at the present time, we have considered the case in which the non- 
minimal gravitational coupling of the electromagnetic field I(T) changes in time only during 
inflation, whereas it does not evolve any more, i.e., /(t/r) = I(t]o), after the instantaneous 
reheating stage following inflation. 
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Finally, we compare our results with the analysis executed in Ref. 55j. The parameter 



of P in our study corresponds to that of v in Ref. 



551 ] . In particular, (3 = 4.2 in the present 



work correlate with v = 2.6. IN 
is obtained for v = 5/2 in Ref. 



ote that the scale-invariant spectrum of the magnetic fields 



55(, In comparison with the analysis in Ref. |55|, for v = 2.6 
the resultant strength of the magnetic field is estimated as 5.4 x 10~ 9 G, whereas in the 
scale-invariant limit the magnetic field would be 1.4410 x 10~ n G. These figures may change 
depending on the assumptions on the reheating stage. Hence, we suppose the sudden (i.e., 
spontaneous) reheating where all the energy density of the inflaton can safely be assumed 
to be released into the energy density of the radiation. In this sense, for the case where 
H R = 10 10 GeV (T R = 8.6 x 10 13 GeV) with the above values such as = 4.2 (i.e., p = 10 
and n = —18.9), we find 



\B(r} , L = 1 Mpc)| = 1.6 x 10" id G . 



(4.12) 



Clearly, this strength can satisfy the scale-invariant limit, namely, less than 1.4410 x 10 G. 



V. CONCLUSIONS 



We have investigated the generation of large-scale magnetic fields in inflationary cos- 
mology in the context of teleparallelism. We have examined a non-minimal gravitational 
coupling of the torsion scalar to the electromagnetic field, which breaks its conformal in- 
variance and hence, the quantum fluctuations of the electromagnetic field can be produced 
during inflation. It has explicitly illustrated that if the form of the coupling is a power-law 
type, the magnetic field with its strength of ~ 10 _9 G and the coherence scale of IMpc at 
the present time can be generated. This field strength is enough to account for the large- 
scale magnetic fields observed in clusters of galaxies only through the adiabatic compression 
during the construction of the large scale structure of the universe without the dynamo 
amplification mechanism. 

Finally, we remark that the resultant field strength of ~ 10 _9 G on IMpc scale is com- 
patible with the upper limit of ~ 2-6 x 10~ 9 G obtained from the observation of CMB 



radiation 



60 



6JJ as well as that of being smaller than 4.8 x 10 9 G from CMB radiation 



on the present strength with scales larger than the present horizon 62|. There also exist 
constraints on the strength of the large-scale magnetic fields from the matter density fluctu- 
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ation parameter as [63 



the fifth science (S5) run of laser interferometer gravitational-wave 



observatory (LIGO) 64J, Chandra X-ray galaxy cluster survey and Sunyaev-Zel'divich (S-Z) 



survey 



651 ]. which are compatible with or weaker than those from CMB. Incidentally, generic 



features of the spectrum of the 
have been investigated in Ref. 



arge-scale magnetic fields generated at the inflationary stage 



661 ] . Moreover, it is also known that from the Big Bang Nu- 



cleosynthesis (BBN), there are limits on the primordial magnetic fields The constraint on 
the current strength of the magnetic fields on the BBN horizon scale ~ 9.8 x 10 _5 /i _1 Mpc, 



where h = 0.7 57J, is smaller than 10 6 G 67|. Furthermore, it is meaningful to note that 
the large-scale magnetic fields with the strength ~4x 10~ n -10~ 10 G at the present time 



can be observed [68| by various future po 



PLANCK [69j,|70j, QUIET 



71 



72j, B-Pol 



arization experiments on CMB radiation, e.g., 
73 ] and LiteBIRD {74}]. If such large-scale mag- 



netic fields in void regions and/or inter galactic medium are detected, the possibility that 
those origin is the quantum fluctuations of the electromagnetic field generated at the infla- 
tionary stage would become higher. Thus, physics in the early universe including inflation 
may be understood through the future detection of the large-scale magnetic fields. 
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